FERMIL AB-Pub-05 /380-T-rev 
BUHEP-05-18 



New states above charm threshold 



o 
o 

>>: 



> 

On 



in 
o 

i 

Oh! 
CD 
43 



X 



Estia J. Eichten/'Q Kenneth Lane, 2 -0 and Chris Quigg 1 '!] 

1 Theoretical Physics Department 
Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, IL 60510 
2 Department of Physics, Boston University 
590 Commonwealth Avenue, Boston, MA 02215 

We revise and extend expectations for the properties of charmonium states that lie above charm 
threshold, in light of new experimental information. We refine the Cornell coupled- channel model 
for the coupling of cc levels to two-meson states, defining resonance masses and widths by pole 
positions in the complex energy plane, and suggest new targets for experiment. 
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I. INTRODUCTION 

In the short time since the Belle Collaboration reported 
sighting the "missing" charmonium state 77^(2 1 So) in ex- 
clusive B — > KKsK~tt + decays pj, new states asso- 
ciated with charmonium have appeared in great profu- 
sion. Apart from the long-sought hd^Pi), recently ob- 
served by the CLEO Collaboration in the decay ^(2S) — > 
■K°h c — > (77)(7??c) 00 all the newly observed states lie 
near or above charm threshold. The identification of lev- 
els above the DD flavor threshold has renewed the impor- 
tance of theoretical studies that couple the cc spectrum 
to charmed-meson-anticharmed-meson channels. 

Motivated by the Belle Collaboration's discovery Q 
of the narrow state AT (3872) — ► ir + ir~J/i/j, confirmed in 
short order by the CDF [|, D0 @, and BaBar Exper- 
iments, we explored the influence of open-charm chan- 
nels on charmonium properties, and profiled the 1 3 D_2, 
1 3 D 3 , and 2^1 charmonium candidates for A"(3872) 
We reaffirmed the single-channel potential-model expec- 
tation 0,^3 that the favored 1 3 D2 and 1 3 D3 candidates 
both should have prominent radiative decays, and noted 
that the 1 3 D 2 might be visible in the D°D*° channel, 
while the dominant decay of the 1 3 D3 state should be 
into DD. We proposed that additional discrete charmo- 
nium levels were ripe for discovery as narrow resonances 
of charmed and anticharmed mesons. Barnes, Godfrey, 
and Swanson have reached similar conclusions using a 
different model for the coupling to open charm 

Experiments have enriched our knowledge of AT (3872), 
but its precise nature remains open to debate. On current 
evidence, AT(3872) is likely to be a J PC = 1++ state. If 
so, the surviving 2 3 Pi charmonium candidate is problem- 
atic, as we discuss in flIV B 21 Not every state associated 
with charmonium need be identified as a cc state; we will 
mention below what some of the interlopers might be. 
An important goal of predicting the properties of char- 



monium states above flavor threshold is to tag new states 
that do not fit the template, and so might represent new 
spectroscopies. 

Four other new particles require confirmation, for each 
has been seen only in a single experiment so far. Belle 01 
has reported F(3943±ll±13) in the decay B -» KcuJ/tp; 
it is a relatively broad state, with total width T = 87 ± 
22 ± 26 MeV. Belle also reports the state A" (3943 ± 9), 
seen in e + e~ — > J/yb + X 13j- It is observed to decay into 
DD*, but not DD, which suggests an unnatural parity 
assignment. The total width is T < 52 MeV. 

Belle has observed a narrow (r sw 20 MeV) state in 
77 -> DD that they call Z(3929 ± 5 ± 2) 0. The 
production and decay characteristics are consistent with 
a 2 ++ assignment, and this state is a plausible x' C 2 (2 3 P2) 
candidate. The most recent addition to the collection is 
F(4260), a broad (« 88 MeV) 1 — level seen by BaBar in 
e + e~ — f 77r + 7r~ J/ib |l5| . with supporting evidence from 
B -> K'J/ybmr 

In this Article, we fix the masses of the 3 1 So and 
2 3 P 2 charmonium states at the positions of y(3943) and 
Z(3931), in order to sharpen our expectations for the 
properties of other states above charm threshold. We 
also extend our previous calculations to include addi- 
tional states at higher masses. We highlight the search for 
narrow structures in charmed meson-anticharmed meson 
pairs, and remark on analogue states to be expected in 
the bb spectrum. 



II. EXPERIMENTAL STATUS 

Let us briefly summarize what experiment has taught 
us about the new levels associated with charmonium, and 
how the experimental facts square with theoretical expec- 
tations. 



A. Properties of r/ c 

partner of ip{2S) now se___ 
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CLEO [HI, and BaBar [±3 experiments. The world- 
average values of the mass and width are M(r]' c ) = 3638± 
5 MeV and T(rf c ) = 14 ± 7 MeV 0. 

The proximity of 0' and ij' c to charm threshold means 
that communication with open-charm channels will al- 
ter the properties of these two narrow levels. As we 
elaborated in Ref. @, the basic coupled-channel in- 
teraction is spin-independent, but the unequal masses 
and quantum numbers of the 2 3 Si and 2 1 So potential- 
model cigenstates mean that their interaction with dif- 
ferent charmed-meson-anticharmed-meson channels in- 
duces spin-dependent forces that affect the charmonium 
states. These spin-dependent forces give rise to S-D mix- 
ing that contributes to the ■0(3770) electronic width, for 
example, and are a source of additional spin splitting. 

In a single-channel nonrclativistic potential picture, 
the 2S hypcrfine splitting is given by M(0') — M(rj' c ) = 
327ra: s |\I/(0)| 2 /9m 2 . Normalizing to the observed IS hy- 
pcrfine splitting, M(J/ip) — M(r) c ) = 117 MeV, we would 
find M (0') - M(rf c ) = 67 MeV, to be compared with the 
observed 48 ± 5 MeV separation. The 2S induced shifts 
calculated in Ref. || draw tp' and rj c closer by 20.9 MeV, 
substantially improving the agreement between theory 
and experiment. It is tempting to conclude that the ip'-rfc 
splitting reflects the influence of virtual decay channels. 



B. CLEO's ho candidate 

Like the rj' c , the spin-singlet L = 1 level of charmo- 
nium has been an object of desire for three decades. 
The barrier to observation has been the absence of al- 
lowed El or Ml transitions from the J PC = 1 states 
readily formed in e + e _ annihilations. The state 
holds a special interest, beyond the tidiness that would 
come from completing the spectrum of narrow states. 
In the nonrclativistic potential-model descriptions that 
have provided such a reliable guide to quarkonium spec- 
troscopy, the central potential is generally taken to con- 
sist of a Coulomb piece with vector Lorentz structure 
(arising from one-gluon exchange) plus a confining term 
with scalar Lorentz structure. For such a potential, the 
hypcrfine splitting between spin-singlet and spin-triplet 
orbital excitations is very small. A significant deviation 
from the expectation that M^Pi) ps (M(1 3 P, 7 )) could 
offer evidence for an unexpected Lorentz structure, or for 
the importance of effects not included in the nonrelativis- 
tic potential-model description. 

An earlier observation of pp — > h c — > ir°r] r n ear the 
l 3 Pj centroid in Fermilab Experiment 760 [21j is not 
sustained by two additional data runs with 2x and 3x 
the luminosity in the successor experiment E835 |22T,l23j. 
E835 does report a narrow (r < 1 MeV) event excess in 
the pp — > h c — > 7r? r channel at a mass M(h c ) = 3525.8 ± 
0.2 ± 0.2 MeV 22]. The signal bin contains 13 signal 
events on a 3-event background, and the signal strength 
is T(h c -> pp)B(h c -> 7 7/ c ) w 11 eV. 

The CLEO Collaboration has examined 3.08 million 



ip(2$) decays, in search of examples of the isospin- 
violating decay to n h c , followed by ir° — » 77 and 
h c — > 7?7c 0- El] ■ They find 150 ± 40 signal counts for 
the inclusive channel, in which r\ c is not reconstructed, 
and 17.5 ± 4.5 signal counts spread over seven hadronic 
decay channels of T] c , for a combined significance exceed- 
ing 5a. In the inclusive sample, the angular distribution 
of photons in the h c -f7j c cascade follows the 1 + cos 2 9 
pattern characteristic of El emission from a spin-1 state. 
The CLEO experimenters infer the ^Pi mass M(h c ) = 
3524.4±0.6±0.4 MeV, with a product branching fraction 
B(V>(2S) -» TT°h c )B(h c -> 777J = (4.0 ±0.8 ±0.7) x 10~ 4 

The l^i level lies 1.0 ±0.6 ±0.4 MeV below the spin- 
triplet centroid, (M(l 3 Pj)) = 3525.36 ± 0.06 MeV (H, 
consistent with theoretical expectations. 



C. Rare decays of 0(3770) 

The color-multipole expansion 0, HM Ull yields sym- 
metry relations for hadronic transition rates between 
charmonium states, but does not enable calculations of 
rates from first principles. By the Wigner-Eckart theo- 
rem for El-El transitions, all the 1 3 D,/ — > mrJ/ip rates 
should be equal (for degenerate 3 Dj states), so a reliable 
determination of the rate r(0(377O) — > mrJ/tfj) is im- 
portant for anticipating the properties of tp2 and 03 . In 
our survey of i?-meson gateways to missing charmonium 
states @, we adopted the value F(l 3 Di — > mrJ/ip) w 
45 keV. 

The CLEO Collaboration has just reported the 
first high-significance (11.6er) observation of the decay 
0(3770) -> 7T+7T-J/0 H3- Their measured branch- 
ing fraction, £(0(3770) -> tt+ti- Jfip) = (1.89 ± 0.20 ± 
0.20) x 10~ 3 , corresponds to a partial width r(l 3 Di — > 
TT+n-J/i') ps (45 ±9) keV. [The CLEO rate is consistent 
with, but much more precise than, the Beijing Spectrom- 
eter (BES) observation [H, £(^(3770) tt+tt" J/0) = 
(3.4 ± 1.4 ± 0.9) x 10" 3 .] A less significant (3.4<r) ob- 
servation of the decay ?/»(3770) — > 7r°7r°J/0 is consis- 
tent with the expectation that the 7r 7r J/0 rate should 
be half the tt + it~ J/ip rate. Accordingly, we shall take 
F(l 3 Dj — > 7T7rJ/0) = (68 ± 15) keV as an improved es- 
timate of the dipion decay rates of the ip2 and 03. The 
increase from our earlier estimate does not alter the ex- 
pectation that the most prominent decay of the 1 3 D2 
state should be 02 — > JXci- 

CLEO has also searched for rare radiative decays of 
0(3770), reporting the preliminary values T(0(377O) — > 
7Xd) = (78 ± 19) keV and F(^(3770) -» 7 Xc2 ) < 
20 keV |23 ■ These observations are in line with expecta- 
tions (r[l 3 Di -» 7Xd) ~ 59 keV and F(l 3 Di -> 7^2) « 
4 keV |8j), and give us no reason to question the esti- 
mated El transition rates for the other ID levels (see 

Aim . 



3 



D. Properties of X(3782) 

As the best studied of the new cc-associated states, 
A (3872) has been subjected to a broad range of diag- 
nostic tests. Upon discovery, A (3872) seemed a likely — 
though somewhat heavy — candidate for tp 2 (or perhaps 
■03), but the expected radiative transitions to Xc states 
have never been seen. For a state that lies approxi- 
mately 50 MeV above charm threshold, the narrow width 
(r(A(3872) < 2.3 MeV) and the absence of a DD 
signal |3(j suggest unnatural parity P = (— The 
BaBar Collaboration has set limits on the existence of a 
charged partner [3j]]. In 1.96-TeV pp collisions, the pro- 
duction characteristics of A (3872) are similar to those of 
0(2S) 0. CDF has determined the fraction of A(3872) 
arising from B decays as 16.9 ± 4.9 ± 2.0% |3^ . 

The 7T7T mass spectrum favors high dipion masses |33j | , 
suggesting a J/tp p decay that would be incompatible with 
the identification of A(3872) — » J/tp as the strong 

decay of a pure isoscalar state. See, however, the interest- 
ing discussion of alternative interpretations in Ref. |34j| . 
Suzuki has suggested [35| that the p-like behavior of the 
7r + 7r~ mass spectrum in these decays may result from 
a dominant decay to toJ/tp slightly off mass shell and 
the small (isospin breaking) u-p mixing. In this case 
the X(3872) could be an isoscalar state as expected in a 
charmonium interpretation. 

Observing — or limiting — the w ir J/ip decay remains 
an important goal [l(i|. An observed J/tp tt + n~ tt° de- 
cay suggests an appreciable transition rate to J/tp uj [3(| . 
Belle's 4.4-cr observation of the decay A(3872) — ► 
J/07 [3(| determines C = +, opposite to the charge- 
conjugation of the leading charmonium candidates. Fi- 
nally, an analysis of angular distributions (3?| supports 
the assignment J PC = 1 ++ [H, but the mass of A(3872) 
is too low to be gracefully identified with the 2 3 Pi char- 
monium state, especially if Z(3931) (cf. § III E 3|) is to 
be identified as the 2 3 P2 level. The x'ci also lacks an 
isospin-allowcd mrJ/tp decay mode, and acquires a sub- 
stantial DD* width if M(xei) > M(D°) +M(D *°). We 
shall examine this hypothesis further in fllV B 21 

If A (3872) is not a charmonium level, what might it 
be [3^? Three interpretations take the near-coincidence 
of the new state's mass and the D°D*° to be a decisive 
clue: an s-wave cusp at D°D*° threshold E3, a D° - 
D*° "molecule" bound by pion exchange |4ll I42I 113 . 
and a diquark-antidiquark "tetraquark" state [cq] [cq] j4J, 
l45j . According to a different four-quark interpretation, 
AT (3872) is a ccqq state organized by the chromomagnetic 
(color-hyperfine) interaction into a (ce)s-(qq)8 pair |46| . 

What distinctive predictions might allow us to put 
these interpretations to the test? On the threshold en- 
hancement interpretation, we should expect bumps at 
many thresholds, but no radial or orbital excitations. If 
pion exchange is decisive, then there should be no ana- 
logue molecule at D S D* threshold. The tetraquark inter- 
pretation suggests that A(3872) should be split into two 
levels, because [cu] [cu] and [cd] [cd] would be displaced 



by about 7 MeV. A first experimental test by the BaBar 
Collaboration is inconclusive 0|. In 61.2 ± 15.3 events 
that fit the hypothesis B~ -> K~X(3872), they deter- 
mine a mass of 3871.3 ± 0.6 ± 0.1 MeV, whereas 8.3 ±4.5 
B° -> A°A(3872) events yield 3868.6 ± 1.2 ± 0.2 MeV. 
The mass difference, 2.7±1.3±0.2 MeV, doesn't yet dis- 
tinguish between one X and two. If diquarks are useful 
dynamical objects, there should be a sequence of excited 
states as well. On the color-hyperfine ccqq interpretation, 
the only analogous level should be a (bb)s-(qq)s state near 
BB* threshold. 

What if the D°D*° threshold is not the decisive ele- 
ment? Hybrid cc-gluon states might appear in the char- 
monium spectrum |47l |. The mass and 1 ++ quantum 
numbers of X (3872) do not match lattice-QCD expecta- 
tions 0, 0, [53 • The valence gluon in the hybrid char- 
monium wave function leads to speculation that the r/J/ip 
decay mode might be quite prominent. The BaBar exper- 
iment [5l| has found no sign of A(3872) — > f]Jjip. More 
generally, it is plausible that the gluonic degrees of free- 
dom should manifest themselves as "vibrational states" 
in the charmonium spectrum |52l l53j| . While theoretical 
estimates of the spectrum of string-vibration modes do 
not reproduce the characteristics of A (3872), we should 
be alert to the appearance of such states elsewhere. 



E. Evidence for states beyond X(3782) 

The other new states associated with charmonium re- 
quire confirmation and elaboration. For the moment, 
experiments have identified four distinct candidates. 



1. F(3940) 

Exploring the J/ip u mass spectrum in decays B — ► 
K J/tp to, the Belle Collaboration observed an enhance- 
ment of 58 ± 11 events near threshold 01 • Treated as 
an s-wave resonance, this state, named V(3940), has a 
mass M(Y) = 3943 ± 11 ± 13 MeV and a total width 
r = 87 ± 22 ± 26 MeV. The mass of this object lies well 
above DD* threshold, so if V(3940) is a charmonium 
state, it would be expected to decay dominantly into DD 
or DD* . A small J/tp to decay rate would complicate the 
interpretation of the product branching fraction, B(B — > 
A'F(3940))S(V(3940) -> J/tpu) = (7.1±1.3±3.1) x 10~ 5 . 
For comparison, the branching fractions for the decays 
B+ -> K+(cc), with (cc) = r] c , J/tp, XcO, Xa, i>' , are be- 
tween 6 - 10 x 10 -4 [23I . This suggests a branching ratio 
of about 10% for V (3940) -> ujJ/tP if F(3940) is a conven- 
tional charmonium state. In contrast, a cc-gluon hybrid 
state mi ght decay preferentially into J/tp or tp' plus light 
hadronsTHEIl. To settle the identity of y(3940) will 
require confirming that it is a single state and distin- 
guishing it from the X (3940), to which we now turn. 
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2. X(3940) 

The Belle [13, and BaBar Collaborations nav c 
observed copious double charmonium production in e + e~ 
annihilations, with strong signals for J/iprj c , J/tpXcOy an< i 
J/ipVc- Searching in the region above DD threshold, Belle 
has observed a new peak at 3943 ± 6 ± 6 MeV, with a 
width T < 52 MeV at 90% C.L. d. 

The new structure, given the temporary name 
A (3940), is seen to decay into DD*, with a branching 
fraction £(A(3940) -> DD*) = 0.96±g; 4 ^ ± 0.22. The 
absence of any signal in the DD channel (for which the 
branching fraction is < 0.41 at 90% C.L.) suggests that 
A (3940) is an unnatural parity state. Its properties 
roughly match those expected for ?7 C (3S), and we shall 
examine that hypothesis further in flIV Al 



3. Z(3930) 

The Belle Collaboration has examined the DD 
invariant-mass distribution in 77-fusion events in which 
the transverse momentum of the DD pair is small. A 
narrow peak of 64 ± 18 events appears at M(DD) = 
3939 ± 5 ± 2 MeV ■ The state, provisionally named 
Z(3930), has a total width T(Z(3930)) = 29±10±2 MeV, 
and T(Z(3930) -> 7 7)B(Z(3930) -> DD) = 0.18 ± 0.05 ± 
0.03 keV, assuming J = 2. The (helicity) angular dis- 
tribution strongly prefers sin 4 6* , which corresponds to a 
spin-2 resonance, over the flat angular distribution char- 
acteristic of sphi-0. Belle accordingly proposes to identify 
Z(3930) as the 2 3 P 2 x' C 2- The measured properties are in 
reasonable accord with those anticipated in Ref. @; we 
probe the 2 3 P2 assignment further in fllVB II 



4. F(4260) 

In a study of initial-state radiation events, e + e~ — > 
(7)7r + 7r~ J/0, the BaBar Collaboration has observed an 
accumulation of events near 4260 MeV in the ir + Tr~J/tp 
invariant mass distribution The excess of 125 ± 23 
events can be characterized as a sing lc J PC = I" 
resonance with mass 4259 ± 8j^g MeV and total width 
r(V(4260)) = 88 ± 23t^ MeV, but a more complicated 
structure has not been ruled out. It is perhaps signif- 
icant that F(4260) lies near the D* + D*~ threshold at 
4224 MeV. BaBar also reports a 3.1cr signal for F(4260) 
in exclusive — > K ± tt + tt^ J/ip decays ^|J. An ex- 
cess of 128 ± 42 signal events corresponds to a product 
branching fraction B(B- -» #-y(4260))S(Y(4260) -> 
tt+tt-J/^) = (2.0 ±0.7 ±0.2) x 10~ 5 . 

The new state has variously been interpreted as "0(4S), 
a ccg hybrid state [5^, [Eg, E3, and the radial ex- 
citation of a diquark-antidiquark state identified with 
A (3872) |4^. Curiously, the n + TT~J/ip enhancement oc- 
curs in the neighborhood of a local minimum in the cross 



section for electron-positron annihilation into hadrons. 
We shall have more to say about possible interpretations 
in XTvU 

III. STRONG DYNAMICS NEAR THRESHOLD 

Near the threshold for open heavy-flavor pair produc- 
tion, light-quark pairs induce significant nonperturbative 
contributions to the masses, wave functions, and decay 
properties of physical QQ states. Lattice QCD calcu- 
lations, once extended into the flavor-threshold region, 
should provide firm theoretical predictions. At present, 
only a phcnomenological approach can offer a detailed 
description of these effects. 

The effects of light-quark pairs can be described by 
coupling the potential model QQ states to nearby physi- 
cal multibody states. In this threshold picture, the strong 
interactions are broken into sectors defined by the num- 
ber of valence quarks. We decompose the full Hamilto- 
nian TL as 

H = Y,'Hn + H I (1) 

a 

where TL n is the Hamiltonian that governs the (cc) + n- 
light-quark sector and Hj is the interaction that couples 
the various sectors. 

The dynamics of the QQ states (with no valence light 
quarks, q) is described by the interaction Tig. In princi- 
ple, excited QQ states with valence or vibrational gluonic 
degrees of freedom would also be contained in the spec- 
trum of 7io- In practice, a simple nonrelativistic potential 
model is used to determine the properties of the bound 
states in this sector. 

The two-meson sector Qq + qQ is described by the 
Hamiltonian Ti.2- As a first approximation, we shall 
assume H2 to be represented by the low-lying spec- 
trum of two free heavy-light mesons. The physical sit- 
uation is more complex. At large separation between 
the two mesons the interactions are dominated by t- 
channel pion exchanges, when they are allowed. For 
states very near threshold such as AT (3872), such pion 
exchange in attractive channels might have significant 
effects on properties of the physical states [6(j • At some- 
what shorter distances, more complicated interactions 
exist and new bound states might arise, e.g. molecu- 
lar states E3, EE EH E3, Ell Furthermore, at small 
QQ separation the QQ + qq sector may contribute sig- 
nificantly to H.2- In particular, as Swanson has empha- 
sized [64(, the pJ/ip and ajJ/ip thresholds lie in the DD 
threshold region. The main landmarks throughout the 
region we consider are shown in Table [fl 

Because current mastery of nonperturbative quantum 
chromodynamics does not sufiice to derive a realistic de- 
scription of the interactions that link the QQ and Qq+qQ 
sectors, we must resort to a phenomenological Ansatz for 
Tij. Following our earlier work we consider here the 
Cornell coupled-channel (C 3 ) model (6^, |6(| |62j for the 
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TABLE I: Thresholds for decay into open charm and nearby 
hidden-charm thresholds. 



Channel 



Threshold Energy (MeV) 



D°D° 
D + D~ 
D°D*° or D'°D° 

D ± D *T 

uj°J/4> 
D+D7 

D*+D*- 

fJ/i> 

D+Dt- or D*+D; 
a°J/ip 

D*+D*- 
A C A C 



3729.4 

3738.8 

3871.2 

3872.7 

3879.5 

3879.6 

3936.2 

4013.6 

4020.2 

4054.7 

4077 

4080.0 

4081.6 

4116.4 

4223.8 

4569.8 



creation of qq pairs. A brief discussion of various other 
models for TLi is contained in the Quarkonium Work- 
ing Group's CERN Yellow Report [63. The open-charm 
threshold region occupies our interest in this study. How- 
ever, analogous effects are present in the bb states near 
BB threshold and cb states near DB threshold. A de- 
tailed comparison of different heavy-quark systems could 
provide valuable insight into the correct form for the cou- 
pling to light-quark pairs. 

The C 3 formalism generalizes the cc model without in- 
troducing new parameters, writing the interaction Hamil- 
tonian in second-quantized form as 



W/ = f ELi / = Pa(r)V(r - r>)p a (r>) : d 3 r dV , (2) 

where V is the charmonium potential and p a (r) = 
^ijj^ {r)\ a \p{Y) is the color current density, with ip the 
quark field operator and A a the octet of SU(3) matrices. 
To generate the relevant interactions, tj) is expanded in 
creation and annihilation operators (for charm, up, down, 
and strange quarks), but transitions from two mesons to 
three mesons and all transitions that violate the Zweig 
rule are omitted. It is a good approximation to neglect all 
effects of the Coulomb piece of the potential in J5J . This 
simple model for the coupling of charmonium to charmed- 
meson decay channels gives a qualitative understanding 
of the structures observed above threshold while preserv- 
ing the successes of the single-channel cc analysis below 
threshold El El- 



A. Mass Shifts 

In the presence of coupling to two-light-quark decay 
channels, the mass to of the quarkonium state ^> is defined 
by the eigenvalue equation 



TABLE II: Charmonium spectrum, including the influence 
of open-charm channels. All masses are in MeV. The penul- 
timate column holds an estimate of the spin splitting due 
to tensor and spin-orbit forces in a single-channel potential 
model. The last column gives the spin splitting induced by 
communication with open-charm states, for an initially un- 
split multiplet. 



State 


Mass 


Centroid 


Splitting 
(Potential) 


Splitting 
(Induced) 


1% 


2 979.9" 


3 067.6 


-90.5 e 


+2.8 


l 3 Si 


3 096.9° 


+30.2 C 


-0.9 


l 3 Po 


3 415.3" 




-114.9 C 


+5.9 


l 3 Pi 


3510.5" 


3 525.3 C 


— 11. 6 C 


— 2.0 




3 524.4 ; 


+0.6 e 


+0.5 


1 :! P 2 


3 556.2" 




+31.9 e 


-0.3 


2 1 S 


3 638" 


3 674 6 


-50.1 c 


+15.7 


2 3 Si 


3 686.0" 


+16. V 


-5.2 


l 3 Di 


3 769.9" 




-40 


-39.9 


1 3 D 2 
l'Da 


3 830.6 
3 838.0 


(3 815) d 






-2.7 
+4.2 


1 3 D 3 


3 868.3 




+20 


+ 19.0 


2 3 P 


3 881.4 




-90 


+27.9 


2 3 Pi 
2 x Pi 


3 920.5 
3 919.0 


(3 922)** 


-8 



+6.7 
-5.4 


2 3 P 2 


3 931 9 




+25 


-9.6 


3% 
3 3 Si 


3 943' 1 

4 040" 


(4015)' 


-66 c 
+22 c 


-3.1 
+1.0 



(3) 



"Observed mass, from Review of Particle Physics, Ref. [20 
'input to potential determination. 
c Observed l 3 Pj centroid. 
^Computed centroid. 

^Required to reproduce observed masses. 
^Observed mass from CLE O la . 
9 Observed mass from Belle ^M. 
^Observed mass from Belle Il3l . 
'Observed 3S centroid. 



Above the flavor threshold, to is a complex eigenvalue. 

The basic coupled-channel interaction Hj given by 
© is independent of the heavy quark's spin, but the 
hyperfme splittings of D and D* , D s and £>*, induce 
spin-dependent forces that affect the charmonium states. 
These spin-dependent forces give rise to S-D mixing that 
contributes to the electronic widths of 3 Di states and in- 
duces additional spin splitting among the physical states. 

The masses that result from the full coupled-channel 
analysis are shown in the second column of Table [H] 
which revises and extends our previously published re- 
sults The new version presented here includes the 
3S levels and takes account of Belle's evidence f° r 
Z(3930), interpreted as a 2 3 P 2 state (cf. 311 E 3)1 . As 
in our earlier analysis, the parameters of the potential- 
model sector governed by 7io must be readjusted to fit 
the physical masses, u>, to the observed experimental val- 
ues. The centroids of the ID and 2P spin-triplet masses 
are pegged to the observed masses of l 3 Di [-0(3770)] and 



6 



2 3 P 2 [Z(3930)], respectively. The assumed spin split- 
tings in the single-channel potential model are shown in 
the penultimate column and the induced coupled-channel 
spin splittings for initially unsplit multiplets are pre- 
sented in the rightmost column of Table [HI The shifts 
induced in the low-lying IS and IP levels are small. For 
all the other states, coupled-channel effects are noticeable 
and interesting. 

An important consequence of coupling to the open- 
charm threshold is that the ip' receives a downward shift 
through its communication with the nearby DD channel; 
the unnatural parity rj' c does not couple to DD, and so is 
not depressed in the same degree. This effect is implic- 
itly present in the early Cornell papers [HE j but the 
shift of spin-singlet states was not calculated there. The 
first explicit mention — and the first calculation — of the 
unequal effects on the masses of the 2S hyperfine part- 
ners is due to Martin and Richard ■ In the framework 
of the C 3 model, we found HEl (<* Table El) that the 
induced shifts draw tp' and rj' c closer by 20.9 MeV, sub- 
stantially improving the agreement between theory and 
experiment. This suggests that the i/j'-7)' c splitting reflects 
the influence of virtual decay channels. In the case of the 
3S system, both the 3% rfc and the 3 3 Si V(4040) com- 
municate with open decay channels, and the C 3 model 
leads to a modest 9-MeV increase in the interval between 
them, as seen in Table HD 



B. Mixing and Properties of Physical States 

The physical states are not pure potential-model eigen- 
states but include components with two virtual (real, 
above threshold) open-flavor meson states. Separating 
the physical state 4" into QQ (^ ) and two-charmed- 
meson components {^2), the resulting decomposition of 
TL by sector leads to an effective Hamiltonian for the QQ 
sector given by 



TABLE III: Wave function fractions (in percent) in the Vo 
(cc) and V2 (D*-*' D^-*') sectors for near-threshold states. 



H + H\- 



1 



H 2 + it 



Hi 



* = 







(4) 



Solving Eqn. Q in the QQ sector determines the mix- 
ing among the potential-model states and the coupling 
to decay channels. Calculational procedures for the C 3 
model based on the coupling Hamiltonian @ have been 
described in detail elsewhere 0, l66l IuTI ]. Even above 
threshold for strong decays, the coupled- channel Hamil- 
tonian approach J3J allows a definition of the mass and 
width of the resonances in terms of the complex eigenval- 
ues and a decomposition of the state into the cc and DD 
components via the associated eigenvector. Of course, 
the width and mass so defined will differ somewhat from 
the resonance mass and width inferred from the observed 
enhancement in any particular open-charm-meson-pair 
channel. 

The two-meson contributions to the wave functions of 
the low-lying cc states are shown in Table IIIII for the C 3 
model. The overall probability for the physical state to 



— 

State 


— — : — 

Wo fraction 






— — ; 

W2 fraction 










LJD 


L) L) + L) D 


D L) 


l 3 Di 


50.03 


u 


17.76 


4.43 


2.42 






d 


15.83 


4.18 


2.38 






s 


1.33 


0.85 


0.79 


l d D 2 


61.10 


u 


0.00 


15.10 


3.30 






d 


0.00 


13.94 


3.21 






s 


0.00 


2.38 


0.96 


PD2 


63.36 


u 


0.00 


12.33 


4.88 






d 


0.00 


11.38 


4.74 






s 


0.00 


2.00 


1.31 


l d D 3 


50.06 


u 


16.51 


1.59 


7.47 






d 


13.17 


1.54 


7.20 






s 


0.25 


0.43 


1.77 


2 d P 


50.03 


u 


16.48 


0.00 


5.23 






d 


17.74 


0.00 


5.04 






s 


4.41 


0.00 


1.06 


2^1 


50.04 


u 


0.00 


19.40 


4.95 






d 


0.00 


18.73 


4.68 






s 


0.00 


1.43 


0.78 


2 d Pi 


50.06 


u 


0.00 


20.82 


3.10 






d 


0.00 


20.41 


2.99 






s 


0.00 


2.00 


0.63 


2 a P 2 


49.95 


u 


4.47 


12.80 


7.74 






d 


4.32 


11.47 


7.17 






s 


0.24 


0.87 


0.97 


3^0 


50.05 


u 


0.00 


20.78 


4.02 






a 


U.UU 




O.IO 






s 


0.00 


1.01 


0.55 


3 3 Si 


50.04 


u 


1.42 


10.12 


11.86 






d 


1.45 


10.24 


11.24 






s 


0.98 


1.58 


1.08 


2 a Di 


50.02 


u 


3.94 


5.34 


12.27 






d 


3.87 


5.49 


11.89 






s 


2.69 


3.03 


1.46 



be in what we have called the Wo (cc) sector, denoted iJ ce , 
decreases as charm threshold is approached. For states 
above threshold the mixing coefficients become complex, 
as shown in Table HVl Configuration mixing effects mod- 
ify radiative transition rates 0, EI, in addition to 
their influence on S-D mixing and spin splittings we have 
already mentioned. 



C. Zweig- Allowed Strong Decays 

Once the mass of a resonance is given, the C 3 formal- 
ism yields reasonable predictions for the other resonance 
properties. We present in Table our estimates of the 
strong decay rates for all the s-, p-, and d-wave char- 
monium levels that populate the threshold region below 
4200 MeV, together with what is known from experi- 
ment. No new experimental information has come to 
light about L > 3 levels; the predictions given in Ref. || 
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TABLE IV: Charmonium content of states near cc threshold. 
The wave function $ reflects mixing induced through com- 
munication with open-charm channels. Unmixed potential- 
model eigenstates are denoted by |n 2s+1 Lj^. The coefficient 
of the dominant unmixed state is chosen real and positive. 
Physical states are evaluated at masses given in Table UTI 

State Principal (cc) Components 

*(l 3 Di) = 0.10e +0 ' 59 " r |2S) + 0.01e +0 - 81 " r |3S) 
+0.69 ID) + 0.10e +om%7T |2D) 

*(1 3 D 2 ) = 0.77 |1D) -0.10 |2D) -0.02 |3D) 

^(l^a) = 0.79 | ID) -0.10 j2D) -0.02 |3D) 

*(1 3 D 3 ) = 0.70 |1D) -0.09 |2D) +0.02e- 98i ' r |3D) 

*(2 3 P ) = 0.11e-°- 38i7r |lP)+0.70|2P)+0.03e +0 - 56i7r |3P) 

*(2 1 Pi) = 0.18e-°- 19i7r |lP)+0.68|2P)+0.07e +0 - 77i7r |3P} 

*(2 3 Pi) = 0.18e~°- 27i7r |lP)+0.68|2P)+0.07e +0 - 69i ' r |3P} 

*(2 3 P 2 ) = 0.18e-°- 09i7r |lP)+0.66|2P)+0.07e +0 - 85i7r |3P} 

I^So) = 0.02e-°- 06l ' r 1 IS) + 0.20e-°' 18l7r |2S) + 0.68 |3S> 
+0.05e+°' 67 " r |4S) 

*(3 3 Si) = 0.02e- a05i7r | IS) + 0.19e-° :mw |2S) +0.67|3S) 
+0.07e +0 - 5i ' 7V |4S) 

+0.04e +0l59<7r ID) +0Me +0 - 59l7V |2D) 
*(2 3 D!) = 0.02e +osoi7V |2S) + 0me-°- 71i7V |3S) 
+Q.Ue-°- 50iw | ID) + 0.69|2D) 



remain current. 

The l 3 Di state ^"(3770), which lies some 40 MeV 
above charm threshold, calibrates the reasonableness of 
our calculated widths. As we noted in Ref. @, our value 
of 20.1 MeV is in excellent agreement with the world av- 
erage, r(V>(3770)) = 23.6 ± 2.7 MeV |U. 

The resultspresented here differ in two respects from 
those of Ref. y|. First, the x' C 2 an d v" masses have been 
fixed, if only provisionally, by experiment. This results 
in shifts of the masses and properties of all the 2P and 
3S states as shown in Tables |n] and IVII Second, the 
present approach allows a more detailed extraction of the 
composition of the charmonium states above threshold. 
These results are shown in Tables HTT1 and llVl 

Along with the current PDG values for the total widths 
of higher 1~~ cc resonances, we show in Table IVl a re- 
analysis of the existing experimental data by Seth [t3 | . 

The natural-parity 1 3 D3 state can decay into DD, but 
its /-wave decay is suppressed by the centrifugal barrier 
factor. Thus the 1 3 D3 may be discovered as a narrow 
DD resonance up to a mass of about 4000 MeV. 

Barnes, Godfrey, and Swanson ^lj have recently re- 
ported extensive calculations of decay widths of higher 
charmonium states in the framework of the 3 Po model 
for quark-pair production. [Shortco ming s of both the C 3 
and 3 Po models are assessed in Ref. j6g. Detailed com- 
parisons (e.g. Ackleh, Barnes and Swanson [7!|) between 
various light quark pair creation models arc highly desir- 
able.] Unlike the analysis presented here, their calcula- 
tion does not resum the effects of coupling to decay chan- 



TABLE V: Open-charm strong decay modes of the charmo- 
nium states near threshold. The theoretical widths using the 
C 3 model |§ are shown. 



State n 2a+1 Lj Mode Decay Width (MeV) 

Experiment Computed 

^(3770) l 3 Di D°D° 11.8 

D + D~ 8.3 

total 23.6 ± 2. 7 b 20.1 

V>(3868) 1 3 D 3 DD 0.82 

total 0.82 

Xc (3881) 2 3 P DD 61.5 

total 61.5 

/i^(3919) 2 L Pi DD* 59.8 

total 59.8 

Xci(3920) 2 3 Pi DD* 81.0 

total 81.0 

X' c2 (3931) 2 3 P 2 DD 21.5 

DD* 7.1 
total 29±10(stat)±2(sys) d 28.6 

r/" (3943) 3 1 S DD* 49.8 

total < 52 e 49.8 

V<(4040) 3 3 Si DD 0.1 

DD* 33. 

D S D S 8. 

D*D* 33. 

, , , / 52 ± 10 6 1 „, 

total ( 88±5 c | 74. 

^(4159) 2 3 Di DD 3.2 

DD* 6.9 

D*D* 41.9 

DsD s 5.6 

D a D* s 11.0 

total { ™ ±»: } 



"Computed from CLEO bra nchi ng fractions ll.'l . 

h Review of Particle Physics L'll . 

c Reanalvsis by Seth |j2- 

d Belle 111. 

e Belle 



nels. The general scale of resulting decay rates is similar 
to those displayed in Tablc[3 For example, Barnes et al. 
determine r(3 3 Si) = 80 MeV and T(2 3 Di) = 74 MeV. 



D. Radiative Transitions 

As Tables ITTT1 and fTvl show, the physical charmonium 
states are not pure potential-model eigenstates. To com- 
pute the El radiative transition rates, we must take into 
account both the standard (cc) — ► (cc)7 transitions and 
the transitions between (virtual) decay channels in the 
initial and final states. Details of the calculational pro- 
cedure are given in §IV.B of Ref. (67]]. There we also 
illuminated the differences between single-channel and 



TABLE VI: Calculated and observed rates for El radiative 
transitions among charmonium levels. C 3 model rates include 
the influence of open-charm channels. Photon energies, in 
MeV, in parentheses. 



Partial Width (keV) 



1 3 D!(3770) -> 
model 
CLEO [74] 


Xc2 7(208) 
3.9 
< 40 


Xd 7(251) 
59 
75 ±18 


Xco 7(338) 
225 
< 1100 


1 3 D 2 (3831) -> 
model 


Xc2 7(266) 
45 


Xd 7(308) 
212 




1 3 D 3 (3868) -> 
model 


Xc2 7(303) 
286 






2 3 P (3881) -» 
model 


Jhb 7(704) 
39 


0' 7(190) 
17 


l 3 Di 7(110) 
1.6 


o3p floors 

model 


1 3 n, 'vi'i ziyi 
5.2 


i J-^ 2 71"'-' / 
2.6 




2 Pi (3920) — ► 
model 


J/0 7(737) 
15 


7(227) 
75 




2 3 P 2 (3931) -» 
model 


l 3 D l7 (157) 
0.4 


1 3 D 27 (95) 
3.0 


1 3 D 3 7(62) 
3.6 


2 3 P 2 (3931) -> 
model 


J/0 7(775) 
14 


0' 7(272) 
120 




0(4040) -> 
model 


2 3 P 2 7(84) 
18 


2 3 Pi 7(132) 
5.4 


2 3 P 7(156) 
20 


0(4040) -> 
model 


1 3 P 2 7(456) 
12 


l 3 Pi 7(495) 
0.4 


1 3 P 7(577) 
0.03 



coupled-channel expectations for the radiative rates. 

Our expectations for El transition rates among spin- 
triplet levels are shown in Table IVll Again the observed 
properties of 0(3770) confirm the reasonableness of the 
C 3 framework. The radiative decay rates for the 2 3 Pj 
and 3 3 Si levels are calculated at the predicted masses. 
All these rates are small compared to the expected open- 
charm decay rates. 



IV. IMPLICATIONS FOR NEWFOUND STATES 



A. A(3940) as 77 C (3S) 

The observation by Belle [3 of a state X(3940) re- 
coiling against the J/tjj system in continuum production 
at the T(4S) region fits the general behavior expected 
for the 77" charmonium level. This state is not seen to 
decay into DD but does decay into DD*, suggesting that 
it has unnatural parity. Figure shows the complex pole 
position of the 77" in the C 3 model as its bare mass is var- 
ied in 10-MeV steps. At 3943 MeV, we estimate a DD* 
width of 50 MeV, which is quite close to the experimental 
upper bound on the total width, r(A"(3940)) < 52 MeV. 
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> 

CD 




3.88 3.90 3.92 3.94 3.96 

Mass (GeV) 



FIG. 1: Variation of the 77" width with mass in the C model. 
The physical mass of the -X(3940) is denoted by a crossed dot. 



80 - 



> 

CD 



40 - 




3.90 3.95 4.00 4.05 

Mass (GeV) 

FIG. 2: Width of the 3 3 Si state as a function of its mass. The 
physical mass of the 0(4040) is denoted by a crossed dot. 

If 0(4040) is assigned to the 3 3 Si level, it is somewhat 
problematic to identify AT (3940) as its hyperfme partner. 
We show in Figure |2l the behavior of the physical mass 
and width of the triplet 3S state as the bare mass is varied 
in 10-MeV steps. To reproduce the observed masses, we 
would require the spin splitting in the charmonium sec- 
tor to be 88 MeV, which is considerably larger than the 
2S splitting, and larger than expected in naive potential 
models. 

If, for the moment, we take the existence, mass, and 
unnatural parity of A" (3940) as well established, then two 
possible sources of the discrepancy deserve futher inves- 
tigation. 

First, it is difficult to determine the true pole posi- 
tion for the 3 3 Si state from measurements of the step 
in R cS = a(e + e~ — > hadrons)/cr(e + e _ — » /J + /-t _ ) |72j . 
The resonance decay widths are determined from fitting 
measurements of AR in e + e _ annihilation to a model 
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for each resonance including radiative corrections. This 
whole procedure is complicated by its dependence on the 
resonance shape, i.e., the expected non-Breit-Wigner na- 
ture of the partial widths for radially excited resonances. 
It may be more useful to produce a model of AR for 
direct comparison with data. Greater resolving power 
between models is possible if the contribution from each 
individual open heavy flavor final state is separately re- 
ported. 

For the C 3 model, the structure of AR(cc) in the 
threshold region was studied in the original Cornell group 
works [HE EE H3 an d later extended to the AR(bb) in 
the threshold region jj^. The structure of AR(cc) and 
AR(bb) has also been studied in 3 P models of quark pair 
creation j7||. There are also some attempts to compare 
the different models [zllzl- 

Second, thresholds corresponding to an s-wave plus a 
p-wave charmed meson are nearby. Unlike the two— s- 
wave charm meson channels, which are p-wave decays 
for the 3S states, these are s-wave decays and thus have 
stronger threshold effects. In particular the lowest mass 
channel D(J = 0;jf = 1/2~)D(J = 0;jf = 1/2+) cou- 
ples to the S 1 ^! T] 1 ' state but not the 3 3 Si V(4040) state. 
This will increase the induced spin-splitting (see Table ITTl 
column 5) between the states. Thus the observed if)" - rj" 
separation will correspond to a smaller potential-model 
splitting (column 4). 



B. The 2P states 



120 




3.94 3.96 3.1 

Mass (GeV) 

FIG. 3: Mass dependence of T(2 3 P 2 



4.02 




1. Z(3930) as X 'c2 



Belle has also observed the new state Z(3930) in two- 
photon production of DD pairs Q. The 2 3 P 2 and 2 3 P 
levels are natural charmonium candidates; experiment fa- 
vors the J = 2 assignment. In constructing Table ITTl we 
adjusted the 2P centroid to give the mass of the J = 2 
state in agreement with observation. At the observed 
mass, we compute T(2 3 P 2 -> DD) = 21.5 MeV and 
T(2 3 P 2 -> DD*) = 7.1 MeV, in reasonable agreement 
with experiment (see Table |Vj We show in Figure |S] the 
variation of the 2 3 P 2 pole position with bare mass in 
10-MeV steps, using Eqn. to fully account for finite- 
width effects. 



2. X(3872) as ^Pi 

If X(3872) is a 1 ++ state, the unique surviving cc can- 
didate would be the 2 3 Pi level [35|. One problem with 
this interpretation is clear: the expected mass (cf. Ta- 
blelTTJ is 48 MeV above the observed X mass. If Z(3930) 
is indeed the 2 3 P 2 level, it is very unlikely that its J = 1 
partner should lie at 3872 MeV. 

If the mass were not a problem, the 2 3 Pi interpretation 
would still be difficult to sustain. The El radiative tran- 
sition to jJ/ip is greatly suppressed because of a cancel- 



-1.0 -0.5 0.0 0.5 1.0 

M - Threshold (MeV) 

FIG. 4: Mass dependence of T(2 3 Pi) -> DD* . 



lation between the cc and DD* contributions. We obtain 
35 and 21 keV for the transitions to J/ip and ip' respec- 
tively from the cc contributions alone. When the effects 
of the charmed meson virtual channels are included, the 
rates change dramatically to 2.4 and 90 keV, respectively. 
This would make the reported observation of the 7 J/ip 
mode highly unlikely. Of course, within the C 3 model 
the isospin breaking between D°D*° and D + D*~ virtual 
states is relatively small. Even if we assume it to be large, 
so that only the virtual D°D*° component is important, 
our analysis suggests that a search for the radiative decay 
to ip' would be valuable. 

Finally, the decay rate of 2 3 Pi — ► DD* increases very 
rapidly above threshold as shown in Figure^ Since the 
total width of V(3872) is less than 2.3 MeV, the possi- 
bility that the 2 3 Pi level lies more than 22 keV (!) above 
the D°D*° threshold is precluded. When the additional 
partial width r ot hcr of the 2 3 Pi into other modes (ggg, 
r )J/ip 1 etc.) is taken into account, the largest tenable mass 
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is lowered to M(D°) + M{D*°) - r ot hcr- Combining this 
result with the large ratio of r(V(3872) -> D 7r°D°) « 
10r(X(3872) -> ttttJ/V') reported by Belle [l| renders 
the 2 3 Pi interpretation very improbable. 

C. V(4260) and other states 

In initial-state radiation events, BaBar 0] has re- 
cently reported a mrJ/tp resonance, y(4260), with a 
width of approximately 88 McV. This state is hard to 
assign within the cc spectrum. To be directly produced 
from a virtual photon, it must be a J PC = 1~~ S or D 
state. The small contribution to AR discourages the 4S 
interpretation, but leaves open the 2 3 Di identification. 
However the C 3 analysis requires that state at 4.16 GeV 
in order to reproduce the structure of AR. The iden- 
tification of the F(4260) with the 2 3 D! state is ruled 
out in the C 3 model because the partial decay width of 
the 2D state (at 4260 MeV) into charm meson pairs is 
125 MeV. This is not the discovery mode for V(4260) 
but exceeds the reported total width. Therefore, this 
state does not have a conventional charmonium interper- 
tation. One plausible explanation for this state is that it 
is a hybrid state I55H5II. The lattice calculations of Juge 
and collaborators |48tl50l| for the QQ potentials with ex- 
cited gluonic modes give a lowest multiplet (Hi in their 
notation) consisting of a quark spin-singlet 1 state and 
a spin-triplet set of states (0 ^ 1 '", and 2 h ). These 
states would be degenerate in the absence of relativistic 
spin-dependent mass corrections. If this interpretation 
were correct, then analogous narrow states would be ex- 
pected in the bb system. 

V. SUMMARY 

The computed properties of the 2 3 P2 and 3 1 So cc levels 
arc in rough agreement with the Z(3931) and X(3943) 



states observed at Belle. We await more statistics and 
confirmation of these results to make a definitive state- 
ment on the identification of these states. 

It is extremely important to solidify what appears to be 
known about V(3872), the best studied of the new states, 
beginning with its 1 ++ quantum numbers. We emphasize 
the need to confirm the jJ/ip decay mode, which fixes 
the C = + assignment, and to determine the X(3872) — ► 
7r°7r° J/0 branching fraction. 

Experiments should continue to search for additional 
narrow charmonium states in neutral combinations of 
charmed mesons and anticharmed mesons. The most 
likely candidates correspond to the 1 3 D 3 , 1 3 D 2 , and 1 3 F4 
levels @- The analysis we have carried out can be ex- 
tended to the bb system, where it may be possible to see 
discrete threshold-region states as well. For any of these 
new states that are threshold bound states their mass 
is primarily determined by the positions of open flavor 
channels. For example, if the X(3872) arises through ad- 
ditional interaction in the two-meson sector, then in the 
bb system we would expect an analog state not far below 
B*B threshold. 
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